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The effect of a high Reynold's number, pressure-driven flow of a compressible gas on the confor- 
mation of an oligomer tethered to the wall of a square-channel is studied under both ideal solvent 
and poor solvent conditions using a hybrid multiparticle collision dynamics and molecular dynamics 
algorithm. Unlike previous studies, the flow field contains an elongational component in addition 
to a shear component as well as fluid slip near the walls and results in a Schmidt number for the 
polymer beads that is less than unity. In both solvent regimes the oligomer is found to extend in 
the direction of flow. Under the ideal solvent conditions, torsional twisting of the chain and aperi- 
odic cyclical dynamics are observed for the end of the oligomer. Under poor solvent conditions, a 
metastable helix forms in the end of the chain despite the lack of any attractive potential between 
beads in the oligomeric chain. The formation of the helix is postulated to be the result of a solvent 
induced chain collapse that has been confined to a single dimension by a strong flow field. 



I. INTRODUCTION 

The effect of flow on the conformation of surface teth- 
ered polymers is an area of study of increasing perti- 
nence as advances in materials science, nanofabrication 
techniques and chemistry are incorporated into the en- 
gines of industry. Whether tethered polymers are used 
to prevent protein fouling in biomedical devices (such as 
through PEGylation of blood wetting devices), anchor 
receptors in lab-on-a-chip enzyme-linked immunosorbent 
assay (ELISA) [lj, replace frictional contacts with vis- 
cous contacts through weeping lubrication 2J or serve a 
structural purpose within the device (e.g. nylon submi- 
cron filters), understanding how a flow field modifies the 
behaviour of a polymer chain can improve the design and 
functionality of these devices pjj. While many of these 
devices consist of polymer layers, where intermolecular 
polymer interactions are common, isolated tethered poly- 
mer chains, where intermolecular polymer interactions 
are absent, can also be present in these devices. 

The behaviour of isolated tethered polymer chains un- 
der pure shear flow was examined theoretically well be- 
fore an experimental apparatus could be devised to study 
the phenomenon. The earliest model consisted of a Brow- 
nian dynamics simulation of a single bead connected by a 
Hookean spring to surface in a linear shear field [lj. Be- 
yond the intuitive prediction that the polymer would ex- 
tend into the direction of flow and adopt a stretched con- 
figuration, this simple model predicted that instead of as- 
suming a steady conformation in the flow field, the teth- 
ered bead would undergo cyclical dynamics that could be 
described as a four part cycle of entrainment, rotation, 
contraction and extension. During the entrainment step 
the bead is entrained in the flow field and the polymer 
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chain is stretched. During the rotation step, the force on 
the bead is balanced by the force constant of the spring. 
In the case of a polymer chain the effective "spring" is 
created by the entropic loss of the polymer chain as it 
is extended and loses configurational freedom. As the 
bead is now fixed at a set distance from the adhesion 
point, it is rotated towards the wall of the channel. In 
the contraction step, as the bead is rotated towards the 
wall of the channel, it rotates out of the flow field and, 
no longer being entrained in the flow field, it retracts 
back towards the wall. Finally, in the extension step, 
as the bead nears the anchor point, the Brownian dif- 
fusion of the bead carries it away from the wall where 
it becomes entrained in the flow field again. These dy- 
namics were subsequently confirmed experimentally by 
studies of tethered DNA oligomers in aqueous solution 
[5J |5] although the resolution of the data was insufficient 
to perform a quantitative analysis of the dynamics. 

Cyclical dynamics also have been predicted by other 
theoretical models [TlllOj. When these dynamics are 
manifest, there is some question as to whether or not 
they are periodic in nature. In a Brownian dynamics 
simulation of a FENE chain in a linear shear field the 
dynamics appeared to be periodic with a period roughly 
equal to IOtr, where tr is the relaxation time of the 
polymer chain [7]. Using a simplified system of three 
beads bound by Hookean springs to the points of an 
equilateral triangle, the study of Holzer et al [9] found 
a Hopf bifurcation in the position of the beads as the 
shear rate was increased. At low shear rates the polymer 
extended into the direction of flow, at intermediate shear 
rates a limit cycle emerged corresponding to periodic os- 
cillations, which became aperiodic at higher shear rates. 
In contrast, cyclical dynamics were not found in studies 
where a rigid bond model was used instead of a spring 
model [TT] suggesting that bond elasticity may be neces- 
sary for this phenomenon to emerge. More recently the 
study of Zhang et al [10] studied cyclical dynamics us- 
ing three methods: Brownian Dyanamics, Lattice Boltz- 
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mann and mixed direct simulation Monte Carlo (DSMC)- 
Molecular Dynamics (MD) method. They found that the 
dynamics were aperiodic and suggested that the cyclical 
dynamics may be caused by a Poisson process like the 
tumbling of untethered polymers under shear flow |12j . 
If the tumbling process is a Poisson process, then the de- 
lay between tumbling events should follow an exponential 
distribution. In addition to the Brownian dynamics ap- 
proach, the effects of flow on tethered polymer chains 
have also been examined using Monte Carlo simulations 
[T3HE], MD PMEE], DSMC [IH] and Multiparticle Col- 
lision Dynamics (MPCD) [HHH]- 

The previously mentioned studies of the effect of a flow 
field on a tethered polymer only considered a pure shear 
flow. However, there is often slip flow within these de- 
vices and there can be an elongational component to the 
flow if the fluid is compressible. The effect of the com- 
pressibility of the medium is particularly important in 
gas flows. As the contraction step of the cycle assumes 
that the fluid has no velocity near the wall, the presence 
of a slip velocity alone alters the dynamics. 

In this paper, the behavior of a tethered oligomer is ex- 
amined in a pressure-driven flow field of a compressible 
gas under two different solvent conditions; a solvent con- 
dition where the oligomer assumes an equilibrium confor- 
mation in solvent that is equivalent to its vacuum config- 
uration (ideal solvent), and a more compact configuration 
(poor solvent) . Under ideal solvent conditions it is shown 
that the cyclical dynamics are restricted to the end of the 
chain and are aperiodic under pressure-driven flow. Un- 
der poor solvent conditions the oligomer adopts a helical 
configuration instead of undergoing cyclical dynamics. 

To simulate the flow on the oligomer, the dynamics 
were studied using a hybrid MPCDgl [23] and MD al- 
gorithm. A MD algorithm was necessary as the oligomer 
model made use of an explicit potential. The MPCD 
algorithm was chosen as the interparticle distance was 
significant relative to the length scale over which dynam- 
ical variables changed and, as a particle based algorithm, 
it was easy to interface with a molecular dynamics algo- 
rithm. The flow profile was modeled on a pressure driven 
compressible fluid flow within a nanofluidic channel using 
slip boundary conditions. 

The following section of the article, section 2, presents 
the oligomer model that was simulated in this article, the 
MPCD algorithm that was used to simulate the fluid and 
the dimensionless numbers that characterize the physical 
significance of these models. The results for the ideal sol- 
vent conditions and poor solvent conditions are presented 
in sections 3 and 4 respectively. The final section of the 
article, section 5, contains a summary of the results and 
the potential implications for practical applications. 



II. THE SIMULATION MODEL 

The essential difficulty with the simulation of the ef- 
fect of channel flow on the conformation of an oligomer 



is the existence of two disparate time and length scales 
within the system associated with evolution of the flow 
profile within the channel and the configuration of the 
oligomer. Due to the elongational component of the flow 
field, a significant portion of the channel needs to be sim- 
ulated to capture the flow profile. In particular, the en- 
tire cross-section of the channel must be simulated, as 
well as a sufficient portion of the channel both upstream 
and downstream of the oligomer to establish the flow 
field and accommodate inlet and outlet artefacts. The 
flow variables such as pressure, density, temperature and 
velocity also require approximately a nanosecond to at- 
tain a steady state. This requires a continuum approach 
or a mesoscopic model. However, to simulate the effect 
of this field on the oligomer, the oligomer must be rep- 
resented by a model of sufficient complexity to capture 
subtle changes in the native conformation, which requires 
integration steps in the femtosecond range. Thus the 
oligomer requires a MD integration approach. As a result 
of these requirements, a hybrid MD-MPCD algorithm [2"2"] 
was used to approach this simulation. The MPCD algo- 
rithm has the advantage of impressive speed of simula- 
tion while preserving hydrodynamics [24 and, being par- 
ticle based, it is straightforward to adapt the model to 
incorporate the MD subsystem for the oligomeric chains. 
The oligomer model itself is simplified from a complete 
molecular description of a homopolymer. It maintains 
the backbone structure of the chain but removes the side 
groups of the chain, incorporating them into the back- 
bone beads. This facilitates the use of a larger integra- 
tion time step and decreases the number of variables that 
could affect the configuration of the oligomer itself, lead- 
ing to results that have general applicability. 

A. The MPCD-MD Algorithm 

Channel flow for a compressible fluid within a 
nanoscopic channel is complicated by several factors; the 
fluid is expanding due to an isentropic expansion as the 
pressure decreases, the fluid is being sheared by frictional 
interactions with the walls of the channel, and the fluid 
is slipping at the walls of the channels. Due to the com- 
plexity of this flow, the flow profile was simulated us- 
ing a MPCD algorithm, and a MD subsystem was intro- 
duced to model the interactions between the fluid and 
the oligomer. 

In order to streamline the standard molecular dynam- 
ics equations, the MPCD algorithm forgoes the use of an 
explicit potential between solvent particles. The domain 
is subdivided into cells of edge length L = 1 nm. At 
discrete time intervals, r c , particles sharing the same cell 
have their velocity rotated about the average cell veloc- 
ity, v, by an angle <f). I.e. velocity of the i th particle in 
the cell, v* is 

v* = R fe (±4>) (v< - v) + v, 
where is the pre-rotation velocity of the particle and 
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Rfc (±</>) is the matrix for a rotation through an angle 
<j>, about an axis k. The magnitude of the rotation is 
kept constant during the simulation, for the simulations 
presented here (f> — tt/A . In three dimensions there are 
two models for the choice of the direction about which the 
rotation is performed. Either the rotation is performed 
about any randomly chosen direction [24] , or the rotation 
is performed about a randomly chosen Cartesian axis of 
the system [25]. The forms of the transport coefficients 
are somewhat different between the two methods [26T 
|2"9"] . The first model has the advantage of assuring the 
proper rotational symmetry for the stress tensor of the 
fluid, while the second method is slightly more efficient. 
For the purposes of modeling channel flow the second 
model is sufficient, thus it is used here. For this fluid the 
transport coefficients are [2"6H3"0"] : 

pk B Tr c (3pg + 2p c + 1 + y/2 (p 2 c + 2< - 0.5) - 2)) 
[L ~~ 2(3 /O 2_2p c -l- v /2(p2 + 2(-0.5)-2)) 
m(p c - 1 + e-P") (1 -2-°- 5 ) 
+ 18r c L ' 

k B Tr c (9p3 + a t ) L 2 (p c - 1) (l - 2-°- 5 ) 



being modeled. I.e., 
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2m(9p3-a t ) 15r c p2 
k B Tr c ( 3p c 



2m \(p c -l + e-Pc)(l-2-°- 5 ) 



where m is the mass of the fluid particle, p is the shear 
viscosity of the fluid, Dt is the thermal diffusivity, D 
is the self-diffusion coefficient of fluid particles, p c is the 
average number of particles per cell and 

a t = (2 ( Pc - 1) (2-°- 5 - 1) - Pc ) {pi + 8 a5 ( Pc - 1)) 
+9.6 (2-° 5 - 1) ( Pc - 1) . 

The algorithm also randomly shifts cells prior to rotation 
and shifts them back following a rotation by a random 
displacement vector S to eliminate the effect of momen- 
tum correlations from particles being in the same cells 
during subsequent rotations [25] [26], however as the mean 
free path length is greater than L, this procedure is not 
strictly required. 

Although there is no local structure or excluded vol- 
ume interactions in the MPCD model fluid, there is still 
momentum transfer between particles. It is possible to 
introduce some excluded volume interactions by modify- 
ing the algorithm to introduce specular collisions between 
particles in cells where the average cell velocity indicates 
that the particles are travelling towards each other [3"T] . 
A closed MPCD system still conserves mass, momentum 
and energy. It can be shown through an H-theorem that 
an MPCD system evolves to a state of maximum entropy 
and, in the low rarefaction limit, reproduces hydrody- 
namics consistent with the Navier-Stokes equation fM\ . 

For the purposes of introducing a physical connection 
between the simulation data and a real fluid, r c was se- 
lected to be equal to the mean collision time for the gas 
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where r = 0.98 A is the particle radius m = 39.943 AMU 
is the particle mass, T = 298.15 K and Pj„ and P out are 
the inlet and outlet pressure. P ou t was held constant for 
all simulations at 82.3 MPa (corresponding to a particle 
density of 20 particles/cell), while Pi n was varied to gen- 
erate a pressure gradient. The MPCD algorithm used 
here was not integrated using the rotational period r e . 
Instead a time step, At = 0.5 ps that was shorter than 
r c was used for numerical integration. As r c v > L, a 
shorter time step was required for the stable application 
of the inlet and outlet boundary conditions. The bound- 
ary conditions operated by modifying the particle dis- 
tributions on a cell-by-cell basis at each time step, thus 
if the distribution changed significantly between inter- 
vals it would produce artefacts in the flow profile. The 
use of a shorter time step also allowed the velocity to 
be truncated at L/At to prevent particles from travel- 
ing more than one cell length in a single time step. As 
At < j \J y^ti very few particles require truncation and 

having a screening range of particle motion simplified 
the construction of the MD susbsystem during each At 
interval. In the simulations of different pressure differ- 
ences presented here, r c ranged from 8At for equilibrium 
simulations to 4Ai for the highest flow rate simulations. 
Measurements that involve the variance of the distribu- 
tion were corrected for the omission of these few high 
velocity particles by multiplying the measured variance 
by the ratio of the variance of the complete Maxwell- 
Boltzmann distribution to the variance of the truncated 
Maxwell-Boltzmann distribution. 

A leap-frog integration scheme with a time step 5t = 
At/ 1000, was used to simulate the oligomer and its inter- 
actions with the fluid |32j. As the fluid particles do not 
interact directly with each other, the algorithm requires 
the integration of the equations of motion of only parti- 
cles that could interact with the oligomer between t and 
t + At. Due to the velocity truncation, the maximum 
distance a particle can be from the oligomer and still in- 
teract is 2L + r c , where r c = L is the cut off radius of 
the fluid-oligomer potential. Hence, in order to interface 
the oligomer model with the MPCD channel, a cell en- 
velope that contains all the particles that are within at 
least 2L + r c of the oligomer is separated from the system 
and integrated with time step St, creating a hybrid MD- 
MPCD system. The overall technique is similar to that 
employed in Malevanets and Kapral [52], but simplified 
due to the velocity truncation. 

Due to the complexities introduced by the elongational 
component of the flow field, mechanisms for generating a 
pure shear flow field such as sliding brick [33] boundary 
conditions or a systemic force field [34] are not appro- 
priate. Instead, the flow field is created by a source-sink 
boundary condition algorithm that establishes a pressure 
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gradient at the inlet and outlet via the addition and re- 
moval of particles from the inlet and outlet of the system 
[3D]. Particles were added from a Maxwcll-Boltzmann 
distribution at a specific temperature to supply thermal 
control. The mean velocity, used for the calculation of 
the thermal velocity, was calculated using a series solu- 
tion to the Navier-Stokes solution for the isothermal flow 
of an ideal gas in a square channel at the inlet, while a 
spatial and temporal local average was used in the out- 
let. The interactions between the solvent particles and 
the wall were modelled by an adiabatic collision operator, 
where the component of particle velocity that is normal 
to the wall is reflected and the remaining components 
are rotated by an angle that is uniformly drawn from 
the [0, 2ir) interval. The algorithm is presented in full 
detail elsewhere [30 , complete with details on the bound- 
ary conditions and implementation. 

The simulations were performed in a 20L x 20L x 
100L channel of square geometry (Fig. [IJ with the 
oligomer chain anchored to the centre of the channel at 
(— 10L, 0L, 50L), with the z-axis representing the direc- 
tion of flow. Initially the oligomer is absent from the 
channel, which is initialized to a pressure corresponding 
to the inlet pressure. The system is allowed to equili- 
brate for 2000 At, corresponding to 1 ns, and the flow 
profile is generated from a washout distribution, where 
the channel is initially filled to the higher inlet pressure 
and the gradient is allowed to form from the drainage 
of particles through the outlet. This has the advantage 
of more rapid equilibration than the inverse fill-up dis- 
tribution, where the system is equilibrated to the outlet 
pressure and particles are added through the inlet to cre- 
ate the gradient, and avoids the creation of shocks within 
the system. At 2000 At, the oligomer is inserted into the 
channel, with a straight configuration extended along the 
x axis, where each bead is separated from the previous 
by the equilibrium bond length. This is done to remove 
some variation between runs and to observe the relax- 
ation time of the oligomer. Each bead is given an initial 
velocity drawn from a Maxwell-Boltzmann distribution 
with y/k B T/m = 1.2456 x 10" 1 L/ At. Any fluid par- 
ticles that have a non-zero potential with the oligomer 
chain at the time of insertion are deleted. 



B. The Oligomer Model 

A critical factor determining the behavior of the 
oligomeric chain under shear flow is the structure and 
potentials used to model the oligomer. As the MPCD 
solvent used in these simulations lacks short range order, 
a similarly simple oligomeric chain was chosen, where 
each repeat unit of the oligomer is replaced by single 
bead and each bead is connected to the subsequent bead 
through a spring potential. This oligomeric chain is 
representative of an oligomer such as polyethylene or 
poly(tetrafluoroethlyene) where there are no large side 
groups on the chain and the repeat unit for the chain 




FIG. 1: Channel Geometry. For these simulations H=20 L 
and L=100 L. The oligomer is anchored at (-10 L, L, 50 L). 

can be reduced to each backbone carbon (e.g. CHi, 
CF 2 ). Whenever physical data were required, the data 
for polyethylene were used to model the oligomer. Addi- 
tional potentials are defined between oligomer beads and: 
other oligomer beads that do not share a bond, MPCD 
fluid particles and the channel wall. For these simulations 
each oligomer bead had a mass that was 0.3512 m and a 
radius of 0.050 L. The length of the oligomer was held 
constant at 100 beads. The parameters for the potentials 
described below are given in table [TJ 

The oligomer chain bond potential was modeled as a 
finitely extendable non-linear elastic (FENE) spring [35] 
with a modification to allow for an equilibrium bond 
length. I.e. the bond potential was: 

FENE ts I I \ Y i,i+l\ \ 

u i i+1 = Kfene a — |ri )i+ ij — In I '- , , 

where = rj+i — is the separation vector between 

the position of bead i and i + 1, a is the equilibrium 
bond length (thus and 2a are the extensibility limits) 
and Kfene is the spring constant for the bond. The po- 
tential is zero at the equilibrium bond length and infinite 
at a separation of cither or 2a. Hence there is a limit 
on the extension of the oligomer, just as there would be 
for a physical bond. In theory this limit is 2a, though, 
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due to the choice of parameter values, the bond length 
stays within the range of 0.9a to 1.1a under the specified 
simulation conditions. This potential was also used to 
tether the oligomer to the wall at the anchor point. 

The interactions between the oligomer beads and the 
fluid particles were modeled by a scaled Lennard- Jones 
6-12 potential of the form: 



u LJ - E ■ 
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where E m i n is the energy at the minimum of the poten- 
tial well and r m i n is the particle separation of minimum 
energy. While the Lennard- Jones potential is a stan- 
dard model for an intcr-molecular potential, the attrac- 
tive portion of the potential when used as an oligomer- 
solvent potential could act to cross-link the oligomer 
when there are no explicit solvent-solvent potentials 
present. Solvent-solvent interactions prevent multiple 
solvent particles from occupying the same volume within 
interaction range of the oligomer. If the attractive poten- 
tial is strong enough to bind a solvent particle, it can bind 
multiple solvent particles to the same site and these sol- 
vent particles can act to cross-link the oligomeric chain. 
Instead of a model that is closer to the physical reality 
of the situation, the model creates a magnified primary 
solvation shell with no secondary structure. Within the 
parameter range chosen to model this interaction, this 
phenomena did not occur as there was no solvent density 
increase in the vicinity of the oligomer. 

The Lennard- Jones potential was used to simulate 
both poor solvent and ideal solvent conditions. A ideal 
solvent is one in which the enthalpy of mixing for the 
solvent and the polymer is zero, and the polymer adopts 
the same conformation that it would in the absence of a 
solvent. In theory an ideal solvent for a polymer would 
be a solvent. Theoretically, under solvent conditions 
the radius of gyration of a tethered polymer is expected 
to have a similar magnitude to an untcthcrcd polymer 
and scale with the length as R gyr — ay/ N/6, where N is 
the number of Kuhn segments in the chain and a is the 
length of the segment |36|. However, unlike a true poly- 
mer, the excluded volume interactions in the core of a 
short-chained oligomer prevent the oligomer from adopt- 
ing Gaussian chain characteristics |37| and the oligomer 
assumes a configuration at least as large as a self-avoiding 
chain [571 l3"B"] . R gyr oc jV 59 , where the proportionality 
constant is specific to the oligomer being examined. As 
such we use the term ideal solvent instead of O solvent 
and state that if the oligomer adopts a smaller confor- 
mation than its vacuum configuration, the solvent is said 
to be a poor solvent, while if it adopts a larger con- 
formation in the solvent than in vacuum, the solvent is 
said to be a good solvent. Hence the solvent conditions 
were estimated by comparing the radius of gyration of 
the oligomer in the solvent, to the radius of gyration of 
the oligomer in the absence of a solvent under equilib- 
rium conditions, as opposed to comparing the size of the 



oligomer to Gaussian chain statistics. The radius of gy- 
ration of an oligomer is written as, 



Rgyr 
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where N is the number of beads in the oligomeric chain, 
q; is the position vector of the i th bead in the polymer 
chain and R gyr ,x, Rgyr,y and R gyr ,z are the x, y and z 
components of the radius of gyration. These components 
are analogous to the radius of gyration and are formed 
using the x, y or z component of q^, i.e. 



R„ 



\ 




In the absence of a solvent, a centimer was found to have 
a radius of gyration of 1.85 ± 0.32 L, while the radius of 
gyration was 1.79 ±0.26 L under ideal solvent conditions 
and 1.21 ± 0.20 L under poor solvent conditions. 

The interactions between beads that were not bonded 
within the oligomeric chain were determined by a trun- 
cated harmonic potential: 

,h . _ / K h r 2 c /2 (l-|r„ |/r c ) 2 |r y | < r c 
\rij \ > r e , 

where is the spring constant of the potential and r c is 
the cut off radius of the potential. Without excluded vol- 
ume interactions, an attractive potential tends to cause 
the chain to favor a collapsed configuration. This con- 
straint coupled with the need for a clean cut-off potential 
led to its selection. 

The interactions between the walls and the oligomer 
were modeled using an exponential repulsive potential: 



-/3d 



d/3ae' 



< d < d c 



where d is distance between the particle and the nearest 
point on the wall, a is the maximum potential between 
bead i and the wall, /3 is the decay constant of the po- 
tential and d c is the cut off radius of the potential. This 
potential was chosen to prevent the oligomer from ad- 
sorbing to the wall of the channel. 



C. Dimensionless Numbers 

Dimensionless numbers are often used in fluid mechan- 
ics and polymer physics to characterize the solvent and 
determine the dominant physical processes within the 
system for the purpose of modeling. Thus far this section 
has presented the simulation models and physical param- 
eters, this section will discuss the physical significance of 
these choices. 
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Parameter 


Value 


Kfene 


43.7216 mL/At 2 


a 


0.154 L 


Emin (Ideal Solvent) 


-3.0483 x 10" 3 mL 2 /At 2 


E m in (Poor Solvent) 


-1.0953 x 10~ 2 mL 2 /At 2 


Tmin 


0.2073 L 


K h 


25.0268 mL 2 /At 2 


r c 


0.308 L 


a 


1.315 x 10 2 mL 2 /At 2 


P 


0.3 L' 1 


d c 


1 L 



TABLE I: Potential parameter values for simulations. 



The fluid velocity profile is governed by four di- 
mensionless numbers, the Knudsen number, K ni the 
Reynolds number, R e , the Mach number, M a and the 
Prandtl number P r . The first three of these numbers are 
related to each other through the equation, 



2 R' 



where C is the ratio of specific heats. The Knudsen num- 
ber is defined as K n = A/A where A = p is the 

characteristic length scale of the system[39]. In essence 
it indicates whether the inter-particle distance is signifi- 
cant over the range that hydrodynamic parameters vary 
significantly. If K n < 0.001 then a continuum method 
such as the Navier-Stokes equation can be used to model 
the flow profile, while a large K n (K n > 10) indicates 
that a molecular dynamics approach is more appropriate 
to simulate the dynamics of the system. In the simu- 
lations with flow presented here, the Knudsen number 
ranged from (8 ± 3) x 10" 4 to (2.08 ± 0.03) x 10~ 2 in 
spatial regions in the channel that are far from the poly- 
mer. In the vicinity of the polymer, local thermodynamic 
and transport properties are not well-defined due to rapid 
changes in the microscopic definition of such quantities 
over molecular length scales. The Mach number is the 
ratio of the speed of the fluid to the speed of sound in 
the fluid, and is defined for a gas as M a = v^J^rp. A 
low Mach number indicates that compressibility effects 
of the fluid can be ignored. For the present simulations, 
the Mach number ranged from to 0.2413 ±0.0006, indi- 
cating that the compressibility effects of the fluid cannot 
be universally ignored. The Reynolds number is the ra- 
tio of inertial forces to viscous forces, and for a square 
channel it is defined as R e = mpvH / ' ji. At high Reynolds 
numbers (R e > 2300 for pipe flow[4"0]). turbulence begins 
to emerge in the flow. Over the course of the simulations, 
the Reynolds number ranged from to 320±1, indicating 
laminar flow. Finally, the Prandtl number represents the 
ratio of momentum diffusivity to thermal diffusivity, and 
is defined as P r = p/(mpDT)- In the simulations ana- 



lyzed here, the Prandtl number had a value of 4.4 ± 0.1, 
indicating that convective mass transfer was dominant 
over conductive mass transfer. 

In the presence of a grafted polymer, another dimen- 
sionless number becomes relevant, the Schmidt number. 
This number indicates the relative strength of convective 
mass transfer processes to diffusive mass transfer pro- 
cesses and also indicates whether hydrodynamic interac- 
tions are gas-like or liquid-like. The Schmidt number, 
defined as S c = p/(pmD), was computed to have a value 
of 0.313 ± 0.003, indicating a gas-like solvent. As the 
flow profile includes the pressure drop in the fluid and 
the SRD fluid described here has an ideal gas equation 
of state, a gas-like fluid should be simulated instead of a 
liquid-like fluid where the compressibility is negligible. 

The diffusion coefficient for the center of mass of the 
polymer, D p , can be written as [41] . 
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B c is the average number of polymer beads per cell and 
nib is the mass of a polymer bead. The ratio of the 
diffusion coefficient of the fluid to the diffusion coeffi- 
cient of the polymer , D/D p = 34 indicating that the 
center of mass of an unattached polymer diffuses signif- 
icantly more slowly than the fluid corresponding to the 
conditions relevant for Brownian dynamics of polymer as 
a whole. Similarly, the ratio of the self-diffusion coeffi- 
cients of the fluid and polymer bead particles, D/Db, can 
be estimated[12] to be around 0.38, which suggests that 
the polymer beads themselves move more rapidly than 
the surrounding fluid particles due to their small mass. 
Under these solvent conditions, the dynamics of the poly- 
mer beads is not Brownian-like due to the small size and 
mass of the individual monomers in the oligomer. 



III. IDEAL SOLVENT CONDITIONS 

The initial response of the oligomer from its initial 
stretched configuration to the fluid is dependent on two 
effects, an initial entropic retraction of the chain towards 
the anchor point as the oligomer attempts to assume a 
random coil configuration, and the drag force of the flow 
field pulling the chain in the direction of flow. At low 
flow rates, the entropic retraction is dominant and the 
oligomer adopts a random coil configuration. As the 
flow rate is increased, the oligomer initially assumes a 
parabolic arc that matches the velocity profile of the 
channel. Then, as the tension in the chain increases, 
the oligomer straightens and rotates until it lies recum- 
bent near the wall of the channel. This manifests as a 
general increase in R gyr , more specifically an increase in 
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the flow component, R„ 



and a decrease in the other 



two components, R gyr ,x and R gyr ,y 



Figure [2] shows R gyr 
and its three components as a function of the shear rate 
at the adhesion point of the oligomer to the channel. 
There appears to be a minimum velocity below which 
the oligomer does not extend in the direction of flow, 
and above this velocity the oligomer rapidly extends to 
a value that is approximately 80% of its contour length. 
Because there are no explicit attractive interactions be- 
tween beads in the chain, there is no impetus for adopt- 
ing a specific conformation and the loss of conformational 
entropy from extension should be the primary resistance 
against oligomer elongation. The entropic resistance of 
a material to physical deformation is the basis of the 
entropic theory of rubber elasticity. Extrapolating this 
theory to a single chain, the relationship between the 
shear rate, 7, and the elongation ratio of the oligomer, 
A, should be l44l. 
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As the A -2 term vanishes very rapidly for large elonga- 
tion ratios, one expects essentially a linear relationship 
between the elongation ratio and the shear rate. How- 
ever the shear force is not the only force acting on the 
oligomer, as the flow field has an additional contribution 
from the slip velocity. At low flow rates, a greater than 
linear increase in the flow rate is expected. At higher 
flow rates, the oligomer reaches its elastic limit as bond 
stretching forces become dominant. Thus a sigmoidal re- 
lationship between the radius of gyration and the shear 
rate is expected and confirmed by the data in Fig. [2j 

Visualizations of the oligomer trajectory revealed that 
there appeared to be two time scales of motion for the 
oligomeric chain; a fast primarily torsional mode of mo- 
tion along the length of the chain that produces little ef- 
fect on the radius of gyration and a slower bulk movement 
of the chain that creates larger fluctuations in the radius 
of gyration. The faster mode manifests as finer noise in 
the data and appears to be periodic, as can be seen in 
the cross-correlation function of velocity components in 
the direction normal to the direction of flow computed 
via the equation 



(V X Vy(t)) = 



1 



(t s - t) <r Va: a v 



(v x (r) - v x ) (v y (t + r) - v y ) dr, 



where t s is the length of the simulation, and Vy~ are 
the time averaged values of velocity components v x and 
v y , and <j Vj . and <r Vy are the standard deviations of v x 
and v y . Periodic oscillations within the x — y plane will 
be evident as a sinusoidal wave within this function. Fig- 
ure [3] illustrates these data as a function of time for two 
data sets with a low flow rate and a high flow rate. The 
oscillations do not seem to be dependent on flow rate, 
nor are they consistent across simulations. They appear 
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FIG. 2: Rgyr and its three components vs. Wall Shear Rate 
at the point of adhesion (0, 10 L, 50 L). 



to be composed of rotations in the frequency range of 
1.75xl0~ 2 to3.25xlCT 2 cycles/A^. A window-smoothed 
spectral analysis of the cross-correlation function of the 
transverse components of the oligomer velocity is shown 
in Fig. [4^i. The spectral analysis of this correlation func- 
tion has a large amount of noise, though the bulk of the 
function can be found in three frequency ranges: the 
to 10 _1 cycles/At range (typically around 1.75 x 10~ 2 to 
3.25 x 10 -2 cycles/At or twice that value), a wide peak 
or cluster of peaks around 1.15 cycles/At and another 
peak around 4 cycles/ At. These higher frequencies cor- 
respond to the normal modes of the oligomeric chain in 
the absence of fluid interactions. They also vanish when 
the velocity of a segment of chain is used instead of a 
single fluid particle (Fig. [ZJd) . Hence these modes do not 
seem to contribute to this torsional reptation. 

The slower mode is aperiodic and uncorrelated, and is 
likely created by local fluctuations in the fluid velocity. 
As the shear rate increases, the motion is restricted to the 
end of the chain. The likely cause of this restriction is 
the increased configurational freedom of the untethered 
end of the oligomeric chain relative to the tethered end. 
As noted below, the deviations in the flow field created 
by the oligomeric chain persist to the end of the chain, so 
it is possible that this effect is assisted by the shielding 
effect of the oligomer wake. 

In the introduction, it was noted that cyclical dynam- 
ics are often observed for simulations of tethered poly- 
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FIG. 3: Cross-correlation function between the x and y ve- 
locity components for the last 50 beads of the centimer. The 
solid line corresponds to a flow rate of (3.26 ±0.10) x 10 4 
particles/At, which corresponds to a wall shear rate of 
(2.3 ±0.8) x 10 -3 A^ 1 and slip velocity of (2.9 ± 0.8) x 
10 -3 L/At. The dashed line corresponds to a flow rate of 
(1.20 ± 0.01) x 10 5 particles/At, which corresponds to a wall 
shear rate of (8.3 ±1.4) x 10~ 3 At -1 and slip velocity of 
(7.6 ± 1.2) L/At. The minimum number of samples used to 
calculate the correlation function was 40000 for the measure- 
ment at 1600 At 

mer chains in shear flow. These were observed in sim- 
ulations with a flow rate greater than (6.0 ± 0.2) x 10 4 
particles/ At, although the cyclical motion only involves 
the latter quarter of the chain, instead of the entire 
chain. These dynamics are illustrated in Fig. [5] It is 
worth noting that these dynamics persist in the presence 
of slip flow, which should hamper the retraction of the 
oligomeric chain towards the anchor point of the chain. 
However, the slip flow is likely the reason why the effect 
is limited to the end of the chain. 

It should be noted that these dynamics do not appear 
to be periodic. As was noted by Delgado et al [7], if the 
dynamics were cyclical this should be apparent in the 
cross-correlation function of the x and z components of 
the radius of gyration. However, this cross-correlation 
function is unremarkable, Fig. [6j It has been suggested 
[TU] , that the tumbling of the end of the chain may be a 
Poisson process. If that is the case, then the distribution 



FIG. 4: Spectral analysis of the cross-correlation function be- 
tween the x and y velocity components. The top graph is the 
function for the 90 th bead in chain counting from the anchor 
bead of the chain. The bottom graph was constructed by av- 
eraging the velocity of the last 50 beads of the chain prior 
to the construction of the cross-correlation function. Both of 
these graphs were calculated from the run with a flow rate of 
(3.26 ± 0.10) x 10 4 particles/ At in[|] 



of the time intervals between tumbling events should fol- 
low an exponential distribution, Fig. [7] Due to thermal 
noise, the tumbling events are difficult to identify except 
by visual inspection of animations of the configurations, 
as shown in Fig. [5j where the coordinated motion of the 
oligomer becomes apparent. The configurations associ- 
ated with the contraction phase of the cyclical dynamics 
are visually distinct and short lived during the simula- 
tions, hence the time intervals between contractions were 
used to define the intervals between tumbling events. At 
lower flow rates the tumbling events are less distinct, 
hence this process was restricted to a few high flow rate 
simulations. The resulting distribution is modeled bet- 
ter by a gamma distribution with a shape parameter of 
2.7±0.7 and a scale parameter of (6.0 ± 2.0) x 10~ 2 A* -1 
than an exponential distribution. As the tumbling of the 
end of the chain requires a certain amount of time and a 
subsequent tumbling event cannot begin until the previ- 
ous event has been completed, there are memory effects 
in the distribution of tumbling events, and the gamma 
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FIG. 5: Illustration of cyclical dynamics in a simulation with 
a flow rate of (1.14 ±0.01) x 10 5 particles/At, which corre- 
sponds to a wall shear rate of (6.7 ± 1.5) x 10~ 3 At' 1 and 
slip velocity of (8.2 ± 1.7) x 10~ 3 L/At. The cyclical dynam- 
ics seem to be confined to the latter quarter of the chain. At 
t = 3411 At the oligomer is extended into the flow field. At 
t = 3457 At the end of the chain has drifted towards the wall 
of the channel. At t = 3477 At the end of the chain has re- 
tracted. At t = 3539 At the end of the chain has drifted back 
into the flow field. At t — 3597 At the flow field has extended 
the chain again. 

distribution is the correct form to observe. 

A relationship between the relaxation time of the 
oligomer and the scale parameter of the gamma distri- 
bution could not be determined. The relaxation time of 
the oligomer, as calculated by fitting an exponential to 
the z-z correlation function for the terminal bead of the 
chain [TO], was very noisy and could best be empirically 
modeled as eexp (-7/C) where e = (4.0 ± 1.2) x 10 2 At, 
C = (1.3 ± 0.9) x 10~ 2 At' 1 and 7 is the wall shear rate of 
the simulation. Due to the quality of the relaxation time 
data, a functional dependence between the relaxation 



time and the exponential parameter could not be deter- 
mined. However, if the tumbling mechanism for the end 
of the chain is similar to the tumbling mechanism for the 
tumbling times of an untethered polymer in shear flow, 
then the distribution can compared to the Rouse model 
of Das and Sabhapandit [12] which models the density 
of tumbling times, Pt{t), as p t (t) — t /a exp (— at/t Q ) 
where t Q is the longest relaxation time of the oligomer. 
For the purposes of comparison, the distribution was 
modeled by an exponential with a parameter value of 
(2.2 ± 0.2) x 10~ 2 At' 1 . Using this model, a « 5 which is 
significantly larger than the value of Das and Sabhapan- 
dit a = 0.325. The value predicted by Das and Sabha- 
pandit is a lower bound on the value of a. Semi-flexible 
polymers are expected to have a higher value of q|45j. 
The cyclical dynamics appear to occur more rapidly than 
the tumbling of an oligomer in free solution because the 
dynamics are restricted to the end of the oligomer instead 
of the entire oligomer, which requires less motion of the 
chain for a cycle to occur, and because the oligomer is 
tethered to the wall of the channel, which increases the 
tension in the chain. 




TimexAt 

FIG. 6: The Cross-Correlation function between the x and 
z components of the Radius of Gyration for the simulation 
conditions in Fig. [5] 

The relative contribution of the slip flow relative to the 
shear flow near the wall can be seen by examining the 
fluid velocity near the wall. By Taylor series expansion, 
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slip contribution is reasonable in the vicinity of the poly- 
mer. It is thus surprising that the contribution of slip 
flow is often neglected in the study of flow on tethered 
polymers, when its presence should be expected in any 
simulation where the Knudsen number is greater than 
0.001. 

Fig. [8] shows a quiver plot of the deviations from the 
normal flow field created by the oligomeric chain. The 
plot was created by comparing the flow profile to a fit of 
the stationary flow profile in the absence of the oligomeric 
chain to the following function: 
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FIG. 7: The distribution of time intervals between cyclical 
dynamics events for a set of simulations with flow rate of 
(1.201 ± 0.007) x 10 5 particles/At, which corresponds to a 
wall shear rate of (7.6 ± 1.1) x 10~ 3 At -1 and slip velocity 
of (8.3 ±1.4) x 10~ 3 L/At. The distribution, shown as a 
bar graph has been fit to an exponential distribution (dashed 
line) and a gamma distribution (solid line). The exponential 
distribution has a parameter of (2.2 ± 0.2) x 10" 2 At' 1 . The 
gamma distribution has a shape parameter of 2.7 ± 0.7 and a 
scale parameter of (6.0 ± 2.0) x 10~ 2 At' 1 . 



the fluid velocity near the abnormal wall that intercepts 
the x-axis at x = —10 L is approximately 



v z (x) sa v s u p + (x+ 10) 



dv z 
dx 



wall 



where v s i ip is the slip velocity and -j^-\ W aii is the shear 
rate of the fluid at the wall. From Maxwell's first or- 
der model for slip flow[46l|47] v s i ip = a\^-\ waU , where 
a = 0.96 ± 0.04 (std. error) is the streamwise mo- 
mentum accommodation and A = 1.11 ± 0.16 L (std. 
error) is the mean free path length. Thus, v z (x) « 
(1.07 + x + 10) ^r\ wa ii- As the oligomer is localized be- 
tween x = —8 and x = —9 L, the slip velocity accounts 
for approximately 33 — 50% of the flow field in which the 
oligomer is entrained. Although the slip flow is not the 
dominant portion of the flow field, it is a significant con- 
tribution. A streamwise momentum accommodation of 
1 is a commonly observed value [57] in channels and this 



(x,y,z) = »< : '( 1 • XX^ ' f*) 



i=0 x j = l 

fc=l l=k 



x 2l y 2k ) 



where a, /3 and 7 are adjustable parameters. The de- 
viations were calculated by subtracting the fitted curve 
from simulation data with the oligomer that had been 
time averaged over 400 At. These deviations were fur- 
ther filtered by removing any deviations that had a mag- 
nitude of less than 4 x 10 -3 L/At, which excluded most 
of the standard noise from the data. However, the fil- 
tering also restricted the analysis to strong flow fields, 
where the threshold value of the filter would not exclude 
the data. The time averaging of the data was neces- 
sary due to the thermal noise (since at any instant the 
thermal fluctuations are on the order of 1.2456 x 10 _1 
L/At). The velocity data were calculated based on par- 
ticle velocities in the MPCD cells. As these cells are L 3 , 
this defines the maximum resolution of the velocity data 
and the hydrodynamic interactions within the system. 
Due to the the choice to simulate pressure-driven flow 
of a gas-like solvent as opposed to a liquid-like solvent, 
the hydrodynamic interactions will be relatively weak. 
Since the oligomer beads are lighter than the fluid par- 
ticles and the thermal fluctuations are relatively large, 
it would require the coordinated motion of several beads 
to produce detectable hydrodynamics within the system. 
As such, a L 3 resolution is quite reasonable. The devi- 
ations show a modification of the flow field around the 
oligomer that begins before the flow contacts the chain. 
The flow deviates in both the x and y directions, like a 
semi-circular dome before the oligomer. There is also a 
decrease in the z component of the velocity. The devia- 
tions persist for the entire length of the oligomeric chain, 
including an acceleration of the fluid after it has passed 
the chain, possibly indicating a low pressure area at the 
end of the chain. For a real fluid, the disruption to flow 
created by a single oligomeric chain would be larger than 
shown in these simulations. While the MPCD collision 
scheme allows for momentum transfer between particles, 
it does not include excluded volume interaction, which 
would be significant as the fluid is compacted against 
the oligomeric chain and has to divert. 
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Flow Deviation Quiver Plot 




FIG. 8: A quiver plot of the deviations in the flow field due to 
the addition of the oligomeric chain for the HF simulation in 
Fig. [1] The data have been time averaged for 400 At and any 
deviations with a magnitude less than 4 x 10 -3 L/At are not 
shown in the plot. The deviations at the end of the channel 
are due to the error in the fit to the flow field which tend to 
be larger at the end of the channel than at the beginning. 





FIG. 9: Equilibrium oligomer configurations, a (left) ideal 
solvent b (right) poor solvent. 



IV. POOR SOLVENT CONDITIONS 

The differences in the dynamical behavior of a tethered 
oligomer are quite stark between poor solvent conditions 
and ideal solvent conditions. Though an ideal solvent 
is often considered a poor solvent, it is the threshold be- 
tween a poor solvent and a good solvent and is considered 
in essence to be a neutral solvent. The equilibrium con- 
figurations between the poor solvent and ideal solvent are 
contrasted in Fig. [9| The key difference in the behavior of 
the oligomer between the poor solvent and ideal solvent 
conditions is that the oligomeric chain tends to adopt an 
extended helical conformation or collapsed globule con- 
figuration under the poor solvent conditions, as opposed 
to a random coil that extends into the direction of flow, 
Fig. [TUj The radius of gyration of the oligomer and its 
three components are plotted versus the wall shear rate 



at the point of adhesion in Fig. 11 and exhibit a similar 
response to the ideal solvent conditions. 

At low flow rates the conformation of the oligomer un- 
dergoes a transition from the extended initial configura- 
tion directly to a collapsed state. As the flow rate in- 
creases, the oligomer typically adopts a metastablc helix 
conformation extending into the direction of flow, which 
collapses to the globule state over time (4000 — 6000 At 
of simulation in some cases). At higher flow rates, the 



oligomeric chain adopts the helix as a stable configu- 
ration for the lifetime of the simulations. The forma- 
tion of metastable helices as a homopolymer collapses 
from a stretched initial configuration has been noted in 
a poor solvent under equilibrium conditions at very low 
temperatures where the thermal fluctuations are small 
enough to allow the metastable state to persist [48], 
but this is first study in which flow is shown to stabi- 
lize this structure. Furthermore, in the study of Sabeur 
the poor solvent conditions were simulated via an attrac- 
tive Lennard- Jones intrapolymer potential that effective 
allows polymer-polymer interactions to stabilize a com- 
pact helical structure, while no explicit attractive intra- 
polymer interactions exist within the helices presented 
here. 

The helix itself exhibits no preference for a left-handed 
or right-handed configuration. There appeared to be 9 
monomers per turn of the helix, which generally grows 
from the end of the chain. Typically there were 3 to 
6 turns in the chain, though that number fluctuated 
over time. The helix itself forms quite consistently, and 
a plot of the cross-correlation function of the x and y 
velocity component of the end of the chain during he- 
lix formation show a cyclical pattern consistent with 
the wrapping of the coil, shown in Fig. [l2j The he- 
lix formed consistently with a twisting frequency in the 
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FIG. 10: Oligomer configurations subject to flow under differ- 
ent solvent conditions. Top: ideal solvent conditions simula- 
tion with a flow rate of (1.20 ± 0.01) x 10 5 particles/At, which 
corresponds to a wall shear rate of (8.6 ± 1.8) x 10~ 3 At" 1 
and slip velocity of 9.1±1.4 L/ At. Bottom: poor solvent sim- 
ulation with a flow rate of (1.56 ±0.01) x 10 particles/At, 
which corresponds to a wall shear rate of (10.1 ± 1.4) x 10~ 3 
At' 1 and slip velocity of (11.8 ± 1.6) x 10~ 3 L/ At. 



(1.9 ± 0.6) x 10 -2 cycles/At range for all the simulations 
where the oligomer adopted this conformation. 

There were no cyclical dynamics observed for the 
oligomeric chain in these simulations. The formation of 
a helix is surprising given the absence of an attractive 
potential between oligomer beads. A visualization of the 
solvent particles in the vicinity of the oligomer did not 
indicate any increase in the solvent concentration within 
the helix or any solvent organization which could have 
been acting as a bridge to stabilize the helix. 

The deviations in the flow field created by the oligomer 
are consistent with the deviations under good solvent 
conditions, Fig. [l3j However, it appears that the largest 
deviations occur at the end of the chain instead of the 
beginning, this is likely due to the helix conformation 
that the chain adopts. It presents a larger cross-section 
against the flow than the chain at the adhesion point and 
hence a larger disruption to the flow profile. 

The helical structures are most likely a result of the 
solvent-induced collapse of the chain from its extended 
initial conformation to a globule being countered by drag 
exerted by the flow field on the oligomer. This hypothe- 
sis is supported by the following observations: The helix 
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FIG. 11: Rgyr and its three components vs. Wall Shear Rate 
at the point of adhesion (-10 L, L, 50 L). 



collapsed to the globule state over time and did not form 
if a collapsed initial configuration was used instead of 
an extended configuration or in the absence of the flow; 
There is no evidence of the solvent stabilizing helical con- 
formations through cross-linking interactions. The chain 
begins in an extended configuration and attempts to col- 
lapse to the globule state. However, the flow field orients 
the chain and exerts a force that keeps the chain ex- 
tended. As the flow field adapts itself to the oligomer, the 
end of the chain is caught in the wake the chain generates 
and it can begin to collapse due to repulsive interactions 
with the solvent. However, it can only collapse within 
the wake of the chain, which is a narrow cone along the 
axis of the extended chain. The torsional reptation of 
the chain, coupled with the linear path along which the 
chain can collapse, causes a local collapse of the end of 
the chain into a helix. In essence, the helix is the re- 
sult of a solvent induced polymer collapse with dynamic 
confinement. 



V. CONCLUSIONS 

The effect of compressible gas flow on an oligomer teth- 
ered to the surface of a nanoscopic channel was examined 
using a hybrid MD-MPCD algorithm. Channel flow was 
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FIG. 12: Cross-correlation function between the x and 
y velocity components for the last 50 beads of the cen- 
timer. The simulation was conducted with a flow rate of 
(1.56 ± 0.01) x 10 5 particles/At, which corresponds to a wall 
shear rate of (10.1 ± 1.4) x 10~ 3 At -1 and slip velocity of 
(11.8 ± 1.6) x 10~ 3 L/At. 



simulated by establishing a pressure gradient at the inlet 
and outlet as opposed to a more artificial method that 
would not replicate both the elongation and shear com- 
ponents of the flow field. While this model ignored the 
effects of explicit fluid-fluid interactions, it did incorpo- 
rate explicit fluid-oligomer interactions and, through the 
MPCD rotation step, momentum transfer between fluid 
particles. As such, it was possible to not only look at 
the effects of the flow field on the oligomer, but also the 
effect of the oligomeric chain on the flow field. Though 
the wake was limited by the gaseous nature of the solvent 
and the lack of excluded volume interactions, it still gave 
an indication of the effect of the oligomer. The simu- 
lations were conducted under two solvent conditions, an 
ideal solvent condition and a poor solvent condition. The 
response of the fluid to the oligomer was consistent un- 
der the two solvent conditions but the response of the 
oligomer to the fluid was quite different. 

The oligomer behaved in a manner consistent with pre- 
vious observations under ideal solvent conditions. The 
oligomer extended into the direction of flow to a value 
roughly equal to 80 % of its extension length. Aperi- 
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FIG. 13: A quiver plot of the deviations in the flow field due 
to the addition of the oligomer chain for the simulation in 
Fig. |12| The data has been time averaged for 400 At and any 
deviations with a magnitude less than 4x 10 -3 L/At are not 
shown in the plot. The deviations at the end of the channel 
are due to the error in the fit to the flow field which tend to 
be larger at the end of the channel than at the beginning 

odic cyclical dynamics were also observed at the higher 
flow rates as the oligomer extended into the flow field. 
The distribution of the intervals between cycles yielded 
a gamma distribution. As the flow has an elongational 
component and a slip component, it is surprising that 
these dynamics persist as it is variation of the shear force 
as a function of distance from the wall that is believed 
to give rise to these dynamics. It is likely that the vari- 
ational component of the flow is large enough relative to 
the constant component of flow to allow these dynamics 
to persist. 

The flow field prompted a transition between a col- 
lapsed globule configuration for the oligomer and an ex- 
tended helix state under poor solvent conditions. Due to 
the initial configuration that extended the oligomer into 
the flow field, there were some realizations of the dynam- 
ics where the oligomer initially transitioned into the helix 
state which then collapsed into the globule state, though 
at higher flow rates the oligomer remained extended. 

These results, in particular under poor solvent condi- 
tions, indicate that the effect of flow rate on the con- 
formation of macromers could be a design consideration 
for lab-on-a-chip devices. In particular, flow induced de- 
naturation of proteins or surface active molecules could 
limit flow rates in some of these devices. However, there 
are also potential advantages to flow induced conforma- 
tional changes. Polymerization of monomers under high 
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shear solvent conditions could result in the creation of 
highly crystalline aligned polymers; it may even be pos- 
sible through flow modulation to create varying layers of 
crystalline and amorphous polymer in the same material 
during the same polymerization. 

The hybrid MPCD-MD algorithm used in this article 
to simulate an oligomer in a nanoscopic channel is only 
one potential application for this simulation technique. 
It would be inefficient to use this method to simulate 
dense polymer networks where the MD subsystem would 
occupy the majority of the simulation space, but the tech- 



nique is practicable for the simulation of channels that 
are sparsely populated with tethered polymers, or to ex- 
amine alternate geometries (e.g. a polymer chain on the 
inside of a right angle bend in a channel) . 
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